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We show that the particular distribution of mass deposited on the surface of a nanomechanical resonator can 
be estimated by tracking the evolution of the device’s resonance frequencies during the process of desorption. 
The technique, which relies on analytical models we have developed for the multimodal response of the system, 
enables mass sensing at much higher levels of accuracy than is typically achieved with a single frequency- 
shift measurement and no rigorous knowledge of the mass profile. We report on a series of demonstration 
experiments, in which the explosive molecule l,3,5-trinitroperhydro-l,3,5-triazine (RDX) is vapor deposited 
along the length of a silicon nitride nanostring to create a dense, random covering of RDX crystallites on the 
surface. In some cases, the deposition is biased to produce distributions with a slight excess or deficit of mass 
at the string midpoint. The added mass is then allowed to sublimate away under vacuum conditions, with 
the device returning to its original state over about 4h (and the resonance frequencies, measured via optical 
interferometry, relaxing back to their pre-mass-deposition values). Our claim is that the detailed time trace of 
observed frequency shifts is rich in information—not only about the quantity of RDX initially deposited but also 
about its spatial arrangement along the nanostring. The data also reveal that sublimation in this case follows a 
nontrivial rate law, consistent with mass loss occurring at the exposed surface area of the RDX crystallites. 


I. INTRODUCTION 

Nanoscale mechanical resonators have proven to be use¬ 
ful tools for chemical detection [1-3] thanks to their incred¬ 
ibly high sensitivity to added mass [4—7] and the ease with 
which their vibrational frequencies can be measured to great 
accuracy. For a device of tens or hundreds of picograms, 
molecules adsorbed onto the surface at the scale of fem- 
tograms or smaller are detectable as shifts in the resonance 
frequencies [8-10]. A serious limitation, however, is that the 
linear relationship between the amount of mass added and 
the change in the observed resonances involves a sometimes 
difficult-to-determine constant of proportionality. This fact, 
which goes unacknowledged in much of the literature, is a ma¬ 
jor impediment to high-accuracy nanomechanical mass sens¬ 
ing. 

For a perfectly uniform distribution of added mass, the con¬ 
stant of proportionality is straightforward to compute: it is a 
simple geometric factor, given by the ratio of the resonator’s 
mode-specific “effective mass” to its true inertial mass [11]. 
Even if the deposition is nonuniform, the situation is still 
manageable so long as the mass distribution is well charac¬ 
terized (and not too concentrated near the nodes of the detec¬ 
tion mode). More typical of a sensing application, though, is 
that the distribution is of arbitrary form and more or less un¬ 
known. In that case, there is no reliable way to extract the 
total adsorbed mass from frequency-shift measurements [12], 
except at the level of an order-of-magnitude estimate. 

To avoid this problem, efforts have been made to concen¬ 
trate mass adsorption to specific sites on a device through 
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complex fabrication [13]. An alternative approach has been 
to employ multimode measurements [14-20], which can pro¬ 
vide some degree of spatial resolution (limited by the propaga¬ 
tion of experimental uncertainties through the “inversion ker¬ 
nel” [21]). The simplest example is that both the size and lo¬ 
cation of a point mass —situated along a resonator with extent 
primarily in one dimension—can be determined from a simul¬ 
taneous measurement of two resonance frequency shifts [22— 
25]. In a slightly different context, multimode measurements 
have been shown to provide single-resonator discrimination 
for the mass added along an array of strongly coupled de¬ 
vices [26], 

The work reported in this paper also takes advantage of a 
multimode framework (described in Sec. Ill), and we find that 
silicon nitride nanostrings [27, 28] are an excellent platform 
for our technique. Indeed, although the original experiment 
by Dohn and collaborators to determine the mass and loca¬ 
tion of a pointlike object was performed on cantilevers [22], it 
was soon recognized that the sinusoidal mode shapes charac¬ 
teristic of nanomechanical strings greatly simplify the analy¬ 
sis [23], As an added benefit, for materials such as stoichio¬ 
metric silicon nitride, the high internal tension that draws a 
doubly clamped beam taut enough to become stringlike leads 
to correspondingly high mechanical quality factors. These 
high values aid in detection [29-3 1 ] and are maintained even 
in the presence of a metallic overlayer added to the device for 
the purpose of functionalization [32, 33]. 

In a previous study on nanostrings, real-time mass sens¬ 
ing was mimicked by carrying out sequential frequency-shift 
measurements of the first two harmonics in conjunction with 
pick-and-place deposition of a single microparticle [23], In 
our work, genuine real-time observations are made: we mea¬ 
sure the first, third, and fifth harmonics simultaneously as a 
function of time while mass sublimates from the device sur¬ 
face. This is possible because we are able to resolve all three 



2 


of these modes using optical interferometry at the midpoint 
of the string with the device entirely under thermomechanical 
actuation. 

We show that measurements of any pair of modes reveal 
not only the instantaneous total mass of molecules adsorbed 
but also their distribution—at least in the regime where the 
distribution is smooth, slowly varying, and roughly symmet¬ 
ric about the resonator midpoint. We provide explicit formulas 
(derived in Sec. Ill and presented in Table II) for the mass co¬ 
efficients of the uniform and nonuniform components of the 
adsorbate distribution in terms of all combinations of the fre¬ 
quency shifts measured in modes 1, 3, and 5—including the 
pseudoinverse result in which three measurements determine 
the two coefficients in the least squares sense—and we show 
that the variance among these estimates can be used as an ef¬ 
fective tool to judge the combined uncertainty due to measure¬ 
ment error and analytical assumptions. 

The multimode measurement technique also tells us about 
the desorption characteristics of the added molecule, here the 
explosive l,3,5-trinitroperhydro-l,3,5-triazine (RDX). In par¬ 
ticular, our analytical model is able to discriminate between 
two situations: (i) where RDX molecules are bound to the res¬ 
onator surface and thus display first-order or Langmuir-like 
rate kinetics; and (ii) where the RDX molecules are bound 
to each other in crystalline formations (born from randomly 
distributed nucleation sites), and sublimation occurs primar¬ 
ily from the crystal interface. As we discuss in Sec. IV, our 
frequency-shift data support the latter interpretation. This is in 
keeping with our expectations, since atomic force microscopy 
studies have already established that RDX forms nanoscale 
crystallites on the silicon nitride surface (see the lower panel 
of Fig. 1 in Ref. 34). 

Our combined experimental and theoretical study of mass 
distribution and sublimation may be useful not only in the de¬ 
sign of future explosive sensors [35, 36] but also in a wide 
range of nanomechanical mass-sensing applications where ac¬ 
curate determination of adsorbed molecular mass is relevant. 



FIG. 1. Optical microscope image showing a 101.7pm nanostring 
with a large number of RDX crystallites physisorbed onto the device 
surface. 


We expect it to have particular importance for single-molecule 
mass spectroscopy [12], 


II. EXPERIMENTAL METHOD 

The nanomechanical devices under study are of a simple, 
doubly clamped beam architecture: 250-nm-thick ribbons, 2- 
3 pm wide and 100-300 pm long. Their out-of-plane vibra¬ 
tional modes oscillate in the MHz range and exhibit quality 
factors exceeding 10 5 . The resonators are fabricated from sto¬ 
ichiometric silicon nitride grown on a sacrificial silicon diox¬ 
ide layer on top of a silicon handle. They are patterned using 
standard optical lithography and reactive ion etching through 
the nitride and then released using a buffered oxide etch. (A 
detailed description of the fabrication process is provided in 
Refs. 26, 31, and 32.) The resulting geometry establishes high 
tensile strain along the length of the beam; as a consequence, 
the device acts as a classical string, with its sinusoidal mode 
shapes [31] enabling the analytical treatments that appear in 
Sec. III. Molecules are adsorbed onto the nanostrings by va¬ 
por deposition of RDX that has been heated and carried to the 
sample chip via a flow of nitrogen gas. It is observed that 
the RDX preferentially adsorbs onto the silicon nitride and 
aggregates in clusters, visible in Fig. 1. We employ various 
masks to gently bias adsorption around the midpoint of the 
string. The weak thermal actuation of the string has no effect 
on where the RDX settles. 

Measurements are performed using optical interferometry. 
Light from a HeNe laser (632.8-nm wavelength) is focused 
onto the nanostrings, with part of the light reflecting from 
the surface of the nanostring and part from the silicon sub¬ 
strate underneath the string. The interference of light there¬ 
fore encodes information on the relative separation between 
the nanostring and the substrate as a modulation of the optical 
intensity. The laser power incident on the string is no more 
than 120 pW, and we have shown elsewhere [31] that there is 
no heating effect. 

We focus the laser light halfway along the nanostring. The 
even-numbered, out-of-plane vibrational modes have a node 
there, and so our measurements are not sensitive to them [11], 
But for the odd-numbered harmonics that we can detect, the 
string center is always the point of maximum amplitude; 
hence we are able to position the laser spot to high accuracy, 
and the signal-to-noise ratio of the measurement is at its local 
maximum. 

The nanostring chip is placed into an optical-access vacuum 
chamber, reducing the viscous damping on the resonators [28] 
and enabling thermomechanical measurements of the nanos¬ 
trings [11], The resulting vacuum (~ 10 4 torr) causes sub¬ 
limation of the RDX from the nanostrings, which is not ob¬ 
served under ambient conditions [35, 36]. The HeNe detec¬ 
tion laser is not resonant with any of the internal levels of 
RDX and therefore its light is not adsorbed by the molecule. 

A high-frequency Zurich lock-in amplifier (model HF2LI, 
capable of demodulating as many as six independent frequen¬ 
cies) is used to isolate up to three harmonics that have large 
thermomechanical displacement at the device midpoint. We 
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track these harmonics as a function of time in order to extract 
useful information about the distribution of mass on the device 
surface and about the characteristics of the RDX sublimation 
process. In practice, the continuous time trace of the optical 
interferometry measurement, z(f), is used to generate a power 
spectral density (PSD), 

S(( 0 \t) = l- ( dxe imi f dx'z(x')z{x' + x), (1) 

1 JO Jt-T 

computed over a sliding time window of duration T ~ 0.85 s. 
The “instantaneous” resonance frequency (O n (t) is obtained 
at each time t by fitting the PSD in the vicinity of its 
nth peak to the usual (nearly-Lorentzian) damped harmonic 
oscillator lineshape, S((0;t) = A n ( 0 „(t )/{[o), 2 (f) — ( 0 2 ]~ + 

\(0 (O n (t) / Q,^\ }, with the quality factor Q,, and the overall 
amplitude A„ also optimized as part of the fit [11, 31]. This 
analysis is simple enough that it can be done concurrently with 
the data acquisition. 


III. FRAMEWORK FOR ANALYSIS 

A. Frequency shift in response to added mass 

We begin by considering the total mechanical energy (ki¬ 
netic plus potential) of a string of length L vibrating in one of 
its normal modes: 


rL 

E n = dxn(x)(olu 2 n {x). (2) 

Jo 

Here, x measures the distance along the string; 0)„ and u„ (x) 
are the angular frequency and displacement profile of mode 
n ; and u (x) is the mass per unit length at position x. A small 
mass perturbation, arising from deposition of a distribution 
of molecules on the surface of the string, leads to a modified 
mass distribution ju(x) —> /i(x) + 5/i(x) and a corresponding 
frequency shift ft),, —> 0)„ + S(O n . To leading order, the varia¬ 
tions 5/i(x) and 3ft),, leave the energy stationary (dE n = 0); 
hence 


_ 2 = / 0 L 3x5/i(x)m 2 (x) 

(On Jq 3x/i(x)h 2 (x ) 

Equation (3) can be understood as a linear relationship 
8 ( 0 n /( 0 „ m/M between the relative frequency shift and the 
ratio of the mass added, m = fgdx 8 fi(x ), to that of the origi¬ 
nal device, M = jgdx fi(x). The constant of proportionality 

_ 1 ( Io dx 8 ji(x)ul(x) \ ( fpdx/i(x) \ 

2 \ fc\dx 8 n{x) j \Jpdxn(x)ul(x)j 

is unique to each mode and depends on how the resonator’s 
own mass and the mass adsorbed on the surface are arranged. 

A variety of other effects can come into play when mass 
impinges on a resonator [37], In devices whose vibrational 
modes are dominated by Young’s modulus terms, a molecular 


covering can induce surface stresses that significantly alter the 
bending stiffness [38] (which, in the case of cantilevers, can 
itself be the basis of a detection scheme [39, 40]). Additional 
mass can also provide a new pathway for energy dissipation, 
so that shifts in the resonant frequencies come about indirectly 
through changes in the mechanical quality factor. Neither of 
these effects is relevant to the devices we are studying. Our 
experiment exploits two properties of nanostrings: (i) they live 
in a high-tension limit where the bending terms are negligible 
and the vibrational modes are almost perfectly harmonic, and 
(ii) they exhibit high mechanical Q that is quite robust to the 
presence of any molecular overlayer. 

For a high-tension string, the mode shape is a sinusoid 
u n (x) = (2/L) 1 / 2 sin()i7z:x/L); and, hence 

^ Sco„ fg dx 8 fi (x) sin 2 (nnx/L) 

(On fg dx jj. (at) sin 2 ( nnx/L ) 

[The reflection symmetry of m 2 imposes the fundamental lim¬ 
itation that 3/1 (x) cannot be distinguished from 3/1 (L — x) 
using frequency-shift measurements alone.] If the unper¬ 
turbed string has a uniform mass distribution /i(x) = M/L , 
then Eq. (5) specializes to 
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( 6 ) 


There are two limits worth emphasizing. In the case of a uni¬ 
form mass deposition profile, with 3/i(x) = m/L leading to 

_ 

(On 2M ’ ^ ; 


the mass added m can be determined from a single frequency- 
shift measurement in any mode. A strongly peaked profile 
represents the extreme opposite case. A point mass m de¬ 
posited at position x m leads to 
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and hence the unknown values of m and x m must be deter¬ 
mined from two frequency-shift measurements in any pair of 
modes: specifically. 
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for modes 1 and 2 (first derived in Ref. 23); and 
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for modes 1 and 3. 
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The general mass-sensing problem, however, is much more 
difficult than either of these two limits and has the character¬ 
istics of an ill-posed inverse problem. The constant of propor¬ 
tionality in the relationship 8 co n /co n m/M is undetermined, 
and there is no way to recover an arbitrary mass profile 8 fi(x), 
except from an infinite number of error-free frequency-shift 
measurements (or, in practice, from knowledge of 8 ( 0 n /( 0 n to 
sufficient accuracy for all n = 1,2,3,... up to a cutoff corre¬ 
sponding to the desired spatial resolution). 

We propose to extract the information from the time evo¬ 
lution of the resonance frequencies in a handful of modes by 
making some reasonable assumptions about 81.1 (x). 


B. Two-parameter mass distribution ansatz 

We suggest a simple, two-parameter form to represent the 
distribution of mass deposited on the string: 


- , mo Km i . f Kx 

S#,w = ir + ir s,n r 


(13) 


Equation (13) is constructed so that mo and m\ have units of 
mass, and the total mass sitting on the string is 

f L ~ . mo Km\ 2L 

/ dx8n(x) = — L+ — -= mo + mi = m. (14) 

Jo L 2L K 

mo is nonnegative, but m\ may be of either sign. m\ >0 de¬ 
scribes a convex mass distribution, whereas —2mo/n < m\ < 
0 describes a concave one. The lower bound on m\, which fol¬ 
lows from 8fJ.(L/ 2 ) ~ mo + Km\ /2 > 0, ensures that the mass 
distribution remains everywhere nonnegative. 

The assumptions behind Eq. (13) are that (i) the distribution 
is symmetric under reflection about the midpoint of the string 
and (ii) the distribution is smooth and slowly varying enough 
that it can be approximated by one component that is uniform 
across the string and another that places mass preferentially 
toward (or away from, when m\ < 0 ) its center. 

Putting Eq. (13) into Eq. (5) gives an expression for the 
frequency shifts in the various modes: 


8 (O n 1 / mo 2 n 2 m\ 

”©7 ~M \Y + 4n 2 -l 


(15) 


Inverting the relationship, we find that the values of mo and m\ 
can be estimated from frequency-shift measurements on any 
pair of modes. For instance, in the case of modes 1 and 2 and 
modes 1 and 3, we find that 
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less useful for distinguishing the components of the mass dis¬ 
tribution, — S(O n /(O n —> (mo + m\)/2M, independent of n, as 
n —> 

A revealing, preliminary application of this style of analysis 
is provided in Table I, where Eq. (17) is enlisted to character¬ 
ize the mass distributed on seven nanostrings of four different 
lengths following various levels of RDX exposure. In each 
instance, the resonant frequencies of modes 1 and 3 are mea¬ 
sured before and after RDX vapor deposition. We find that 
the estimates for m\ exhibit both positive and negative sign 
(meaning convex and concave mass profiles), and in some 
cases show magnitude \mi\ comparable to mo itself. This is 
evidence for a quite substantial variation in how RDX settles 
along the length of the string from one experiment to the next, 
and it confirms that a mass distribution close to uniform is 
only rarely the outcome of our sample preparation. The re¬ 
sults in Table I are fully consistent with our intentional biasing 
of the vapor-deposition process. 

How confident should we be in this analysis? One issue 
that arises is how to quantify the uncertainty arising from our 
choice of ansatz. The coefficients appearing in Eqs. (16) and 
(17) are specific to our choice of the sine function in Eq. (13). 
As a test of the robustness of our assumption, it will be helpful 
to check the results against an alternative functional form. We 
therefore also try a symmetric polynomial 
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It is also important is to look for consistency between the re¬ 
sults achieved with measurements on different sets of modes. 
As discussed in Sec. II, the laser in our optical detection sys¬ 
tem is parked at the string midpoint, so only the odd-numbered 
modes are measured. Table II lists all the possible mass deter¬ 
minations using frequency shifts in modes 1, 3, and 5 (exact 
determinations with the three available pairs and one overde¬ 
termination with the full triplet of modes), under the assump¬ 
tion of a mass profile following either Eq. (13) or Eq. (18). 

A reasonable final value of the fractional mass added comes 
from averaging those eight estimates: 

m niQ + mi 

M ~ M 

(5t»i t )©3 8 co 5 

= 0.108 031-- - 0.534596 —- - 1.573434—-. 

©1 ©3 ©5 

(19) 


The total uncertainty, (A e [m/M]' + A a [m/M] 2 ) 1//2 , is taken 
to be a combination of the experimental errors arising from 
imperfect knowledge of the relative frequency shifts. 
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According to Eq. (15), higher modes (larger n) exhibit de- and the uncertainty due to the choice of ansatz, which we ap- 
creasing differentiation in the m\ coefficient and so become proximate by the spread around the mean estimate given in 
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Length 

(pm) 

Width 

(pm) 

M 

(Pg) 

Trial 

fl 

(MHz) 

h 

(MHz) 

5/i 

(kHz) 

S/ 3 

(kHz) 

mo/M 

m\/M 

m = mo +m 1 
(Pg) 

m u 

(Pg) 

m u /m 

101.07 t 

2.83 

218.10 

1 

2.5823 

7.7928 

-71.083 

-210.511 

0.0506 

0.00340 

11.8 

12.0 

1.02 

172.60 

2.05 

269.80 

1 

1.4871 

4.4745 

-0.100 

-0.237 

9.56x 10~ 6 

93.7 x 10~ 6 

0.0279 

0.0363 

1.30 




2 

1.4885 

4.4785 

-1.708 

-11.103 

0.0139 

-0.00874 

1.41 

0.619 

0.441 




3 

1.4898 

4.4821 

-13.045 

-48.873 

0.0363 

-0.0141 

5.99 

4.73 

0.788 




4 

1.4934 

4.4931 

-2.115 

-7.752 

0.00554 

-0.00203 

0.947 

0.764 

0.807 

216.34 

2.56 

422.30 

1 

1.2104 

3.6304 

-9.891 

-26.318 

0.00827 

0.00605 

6.05 

6.90 

1.14 

309.49 

2.72 

641.88 

1 

0.8345 

2.5067 

-2.163 

-5.641 

0.002 19 

0.00224 

2.85 

3.33 

1.17 


TABLE I. Results of a sequence of RDX depositions onto nanostrings spanning four different geometries, organized by string length and 
deposition trial. M is the total predeposition mass inferred from the material density and device volume. Experimentally determined frequency 
values are reported in columns 5-8, viz., the initial mode frequencies of the first and third harmonics, f\ and /j, and the shifts in those 
frequencies after the mass deposition, 8f\ and 5 / 3 . The ratios mo/M and m\/M are obtained from Eq. (17). Traditionally, one would use 
the shift of the first harmonic alone in conjunction with Eq. (7) to determine the adsorbed mass; we label this quantity m u , since it is based 
on the assumption of a uniform mass distribution. Comparison of m a with the corresponding value m = m 0 +m 1 , coming from our improved 
two-component analysis, reveals that the two estimates can differ quite substantially. Fluctuations in the ratio m a /m indicate that the traditional 
approach routinely over- or underestimates the adsorbed mass in the range of 10% to 50%. This quantitative comparison demonstrates the 
absolute importance of a multimode analysis for accurate mass sensing. The dagger (f) marks the one experiment in which the RDX deposition 
process is completely unbiased; the others involve some degree of masking in order to engineer a nonuniform mass distribution. 


Eq. (19): 
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The notation A f , [i) 0 )„/o)„] in Eq. (20) is meant to indicate the 
error in the relative frequency shift, understood as the ratio of 
two inexact quantities computed according to A e [5o)„/o )„] 2 = 
(A e [8cOn]/co„) 2 + (A e [co„]Sca n /co^) 2 . 


C. Sublimation model 

The preliminary results appearing in Table 1 are obtained 
from two discrete measurements of the devices’ resonance 
frequencies, before and after RDX exposure. We have the ca¬ 
pability, however, to make ongoing measurements of the res¬ 
onance frequencies; these change continuously in time as the 
adsorbate molecules are removed from the nanostring because 
of the vacuum environment. The main thrust of this paper is 
to show how we can exploit this much richer data set. 

Therefore it is useful to sketch out a model of how the mass 
on the string evolves with time in our experimental setup. The 
basic assumption is that the molecules are either residing on 
the string (m) or existing as vapor in the chamber (m v ). RDX 
sublimates from the string surface at a uniform rate a, and 
there is an incoming flux of molecules returning to the surface 
due to intermolecular collisions, denoted by /j. Finally, there 
is y, the rate at which the chamber is being evacuated. This 
picture leads to a coupled pair of rate equations: 


As a consequence, there are two rate constants, given by the 
eigenvalues of the matrix. When y = 0, the eigenvalues are 
Ai = 0 and A 2 = O' I /j; hence, the system equilibrates at a 
rate a + /3 to a steady state with m v /m = a/ /3. On the other 
hand, when y is a fast rate (i.e., larger than both a and /j, 
which is the experimentally relevant situation) we find that 


Ai = 
^2 = 


aB exB(ex-B) 3 . 

- a+ Y + f +Q(y 


a/3 a/3(a-/3) 

y 


^' +0(7 3 ). 


(23) 


The full time dependence of m and m w is then given by 


( m{t) \ = yw(0)+/3w v (0) / y\ {a _ a p /y)t 

\m v (t)J y 2 + a/3 Vaj 

ym v (0)-a;n(0) /-/3\ {y+ p +a p /y)t 
y 2 + a/3 \Y ) ' 

For times t +> (y+ /3 + a (i /y) ~ 1 , the mass on the string de¬ 
cays according to 


m(t) 


m{ 0) + (/3/y)w v (0) a(1 _p /r)t 

l + a/3/y 2 


(25) 


Moreover, if the chamber is being very aggressively evacu¬ 
ated, then the behavior looks like 


m{t)=m{ 0)e~^ ttt , (26) 

with a e ff = a — a/3 /y ss a close to the intrinsic sublimation 
rate of the adsorbate, and we can safely proceed as if m v « 0 
at all times after the pump is activated. 

A refinement of this model is to consider the possibility that 
the rate of mass loss goes as a fractional power of the current 
mass load. This may be appropriate here since the RDX on 
the surface is known to aggregate, rather than arranging in a 
smooth monolayer. Sublimation in that case is likely to be 
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Ansatz 

Modes 

mo/M 

m\/M 

Eq. (13) 
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TABLE II. Estimates of the uniform (in o) and nonuniform (mi) mass components as a function of the relative frequency shifts. A solution to 
the system of equations is given for each pair of modes; also shown is the least-squares, pseudoinverse solution involving all three modes. 


at least partially controlled by loss from the surface area of 
the RDX crystallites (A ~ nr' 3 ), which puts a lower bound 
of 2/3 on the effective scaling exponent. If the distribution of 
crystallites along the device is roughly uniform, then we can 
account for this situation as follows; 

^ = -ocm*(m/m*)^. (27) 

This description requires that we introduce a new material- 
specific mass scale m k and a phenomenological exponent 0. 
We anticipate a value 2/3 < 0 < 1, with conventional expo¬ 
nential decay recovered at the upper end of that range. When 
0 y/ 1, we find that 


dm 


- d(m l ‘h = —ami ^dt. 

1-0 3 ’ 


This leads to 


m(t) = m( 0) 



1/(1 


(28) 


where, t r = [m(0)/ m t ] 1 ^/a(1 — 0) is the finite removal time 
after which all RDX has left the nanostring. 

We now emphasize an important feature that distinguishes 
between the two kinds of behavior. When 0 = 1, the decay is 
exponential at a constant rate 


1 m(t) 

-In —v = a. 

t m( 0) 


(29) 


For a nontrivial value of the exponent, however, the mass loss 
is characterized by an apparent decay rate that increases over 


time: 


Am'-M 


1 


t m( 0) (1 — 0)/ 

1 


In 


l — — 

tr 


t r 


(1 — 0)f r 2f r 3 1 } 

I m* i 1 

= a( W) 


(1 0 )at ( y- 

' 2 {m( 0 )J 


(30) 


Finally, a more complete description must account for a 
mass distribution that varies along the length of the string. 
We treat the deposited mass (5/1 (x,/) as a space- and time- 
dependent field subject to the governing equation 


= Ol) 

Here, /i* is a stand-in for m k jL. In the previously considered 
situations, where either 0 = 1 or the initial mass distribution 
is uniform, the mass distribution simply shrinks away while 
preserving its overall shape. That is not generally true: 


(5/1 (x,/) 


S/tM) 1 ^ 


(l-0)a//* 1 0 / 


i/(i- 0 ) 


(32) 


Equation (32) also makes clear that, unlike in Eq. (28), (5/1 
reaches zero at a removal time that is different at each point 
along the string. 

With Eq. (32) in hand, it is straightforward to obtain the 
time-dependent total mass 

m(t) = [ dxdji(x,t) (33) 

Jo 

or any of the relative frequency shifts 


S(Q„{t) 

(On 


1 f L 

— / dx8u(x,t) sin - 
M.J o ' 


/ llTlX \ 

v~w 


■ 


( 34 ) 
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The latter is simply Eq. (6) with Sp(x) replaced by 8 fj.(x,t). 


IV. RESULTS 


We now apply the models developed in Sec. Ill to our ex¬ 
perimentally acquired frequency shifts. The results reported 
here are taken from measurements on a single device of length 
310 pm and mass 641.88 pg. The device’s first, third, and fifth 
modes are tracked over the course of 4 h. 

Our approach is to match Eq. (34), for the cases of n = 
1,3,5, to our time traces of the resonance frequencies. The 



Time (s) 



FIG. 2. Multimode frequency shifts as a function of elapsed time 
and the resulting RDX mass extraction. In the upper panel, each 
experimental data point represents a relative shift in the resonance 
frequency of one of modes 1, 3, or 5, obtained from the finite-time- 
window PSD [as per Eq. (1)]. (The time-window resolution is fine 
enough that the experimental points bleed together into nearly con¬ 
tinuous curves.) For visual clarity, the mode 3 and mode 5 data 
sets are translated upward by 0.000 25 and 0.0005, respectively. The 
three solid lines are the result of matching our theoretical model [viz., 
Eq. (34) in combination with Eqs. (32) and (35)] to the experimental 
observations. The lower panel shows the corresponding mass. Two 
estimates are presented: one is a discrete set of values calculated 
at each instant according to Eqs. (19-21); the other is a continuous 
curve [produced numerically with Eq. (33)] based on the global fit 
achieved in the upper panel. The two estimates are found to be in 
excellent agreement. 


integrals are carried out numerically with Eq. (32) serving 
as the model for dll (x,t). In the spirit of our original two- 
component ansatz, we choose to express the initial mass dis¬ 
tribution in this slightly more expressive form: 


<5 ju(*,0) = — + —N(p) 


X - 1 - X - 
L\ L 


(35) 


The exponent p adds some additional flexibility with regard to 
the shape of the nonuniform contribution [and permits a rough 
interpolation between Eqs. (13) and (18)]. The normalization 
factor 


N{p) 


2 1+w T(3/2 + p) 

V” r(i + P ) 


is chosen to preserve the property that m = JdxSfi(x, 0) = 
mo + mi. 

The set of independent variational parameters mo, mi, p, 
0, and C=(l — ° is simultaneously optimized us¬ 

ing the nonlinear least-squares Levenberg-Marquardt algo¬ 
rithm [41, 42]. The number of degrees of freedom—only 
five — is remarkably small relative to the large number (« 
3 x 7000) of data points. We find that the total mass on the 
string is mo +mi =2.871 (2) pg and that the decay exponent 
has a value (j) = 0.826(6). The resulting high-quality fit is 
displayed in the upper panel of Fig. 2. The lower panel of 
the same figure shows the mass extracted at each instant from 
frequency-shift measurements according to Eq. (19), with er¬ 
ror bars on the data points obtained from the quadratic mean 
of Eqs. (20) and (21). There are a handful of points in the plot 
where the uncertainty estimate blows up, but this amounts to 
just a few blips in the more than 7000 data points. The solid 
line superimposed on the data points is produced by integrat¬ 
ing [via Eq. (33)] the evolving mass profile that emerges from 
the fitting procedure [Eq. (32) with globally optimized param¬ 
eters]. The consistency between these two approaches gives 
us confidence that our estimate of the total adsorbed mass is 
highly reliable. 

Figure 3 shows snapshots of the mass profile at various 
times in the experiment. Such a determination of the real¬ 
time mass distribution may have important applications, such 
as the study of diffusion of molecules along the surface [43], 
or atomic layer reconstruction [44], of a nanomechanical res¬ 
onator. 

Figure 4 indicates the effective instantaneous decay rate, 
calculated as per Eqs. (29) and (30). One can see clearly that 
the rate is not constant in time—and hence it is incompati¬ 
ble with pure exponential decay. Rather, it seems to increase 
steadily. Moreover, between 8000 s and 11 000 s, it turns up in 
a way that is consistent with the fractional power of 0 found 
in our fit to the frequency-shift measurements. It is worth reit¬ 
erating that this implies that sublimation occurs from the sur¬ 
face of crystallites—and not uniformly from the surface of 
the device—in agreement with optical (see Fig. 1) and atomic 
force (see Ref. 34) microscopy. 

Something further we have checked is whether the fitting 
can be meaningfully improved by allowing a —> a(x) to vary 
along the length of the string. One might imagine, for in¬ 
stance, that the sublimation rate has a uniform contribution 















8 



FIG. 3. Profiles of the mass distributed along the nanostring of 
length 310 |rm and mass 641.88 pg at representative moments in time. 
Shown here are plots of Eq. (32) evaluated with optimized variational 
parameters at the appropriate time t. Accompanying each distribu¬ 
tion is a label giving the time of the snapshot and integrated weight 
under the curve. 



FIG. 4. The instantaneous rate at which RDX desorbs from 
the nanostring clearly increases over time. The data shown are a 
logarithmic rescaling of the data appearing in the lower panel of 
Fig. 2. A linear fit over the time interval [250 s, 7750 s] produces 
a(t) = (Xq + a\t = 0.00032865(6) + 1.445(1) x t. Substantially 
better agreement is achieved with a function having the form of 
Eq. (30) with exponent (/) = 0.8278 (fully consistent with the value 
(j) = 0.826(6) used in the upper panel of Fig. 2) and removal time 
t r = 17247 s. 


controlled by the partial pressure of RDX in the evacuated 
chamber and a locally varying component that comes about 
in some unspecified way from the influence of the mechanical 
motion. We have carefully considered various ad hoc mod¬ 
els, but we find consistently that allowing spatial variation in 
a(x) does not improve the fit enough to justify the additional 
variational degrees of freedom [45], 


V. CONCLUSIONS 

Taking advantage of the sinusoidal mode shape of silicon 
nitride nanostrings under high tension, we develop analyti¬ 
cal models for the multimode frequency shifts that occur as a 
result of mass deposited onto such devices. We apply those 
models to real resonators in the lab and demonstrate our abil¬ 
ity to make reliable estimates of the total amount and distribu¬ 
tion of adsorbed molecules. This work provides experimental 
determination of a nonuniform and nonpointlike mass distri¬ 
bution. 

Our multimode technique produces reliable, time-evolving 
estimates of the total mass adsorbed on the resonator with a 
sensitivity on the order of a few femtograms, even in the ab¬ 
sence of detailed knowledge of the initial mass distribution. 
The time trace of frequency measurements in a handful of 
low-lying normal modes is sufficient to reconstruct that miss¬ 
ing information and thus fix the otherwise unknown constant 
of proportionality [appearing as Eq. (4)]. 

Notably, the estimate of the mass added, as determined by 
our analysis, is often found to disagree with the value arrived 
at by the more traditional approach, which relies on a fre¬ 
quency shift in a single mode and incorporates no understand¬ 
ing of how the added mass is arranged on the device. This 
reveals the extreme importance of such an analysis to accu¬ 
rate mass sensing with nanomechanical resonators. 

In addition, because the adsorbed molecules sublimate from 
the surface of the nanostrings in our experimental system, we 
have had an opportunity to analyze their desorption charac¬ 
teristics. Real-time measurements allow for significantly im¬ 
proved estimates of the mass distributions and reveal the evo¬ 
lution of those distributions over time. This may prove a use¬ 
ful tool in studying molecular diffusion [43] or the formation 
of monolayers onto nanomechanical resonators [44], 

Finally, multimode analysis also leads to insights into the 
form in which the molecules are removed from the nanos¬ 
tring surface by the vacuum. Specifically, we are able to de¬ 
velop sublimation models that account for uniform-rate loss 
from the surface and local-mass-dependent loss from the sur¬ 
face area of a crystallite. Only the latter is consistent with the 
observed frequency shifts and the microscopy of the devices. 
This may have important applications in real-world sensing 
of explosive molecules such as RDX, for example in airport 
passenger and luggage screening. 
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